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I. INTRODUCTION 


The goal of this thesis is to study the combustion 
behavior of porous graphite . The mathematical model devel- 
oped by VATIKIOTIS (1) is a transient, one-dimensional model 
of a porous medium consisting of spherical particles , or 
cylindrical fibers of graphite 

The model covers several aspects : mass transfer, heat 
transfer, combustion , airflow , and temperature dependency 
of thermophysical properties. Geometric parameters of the 
porous medium include porosity (p) and medium thickness (L). 
The porous medium has interconnected pores which permit air 
flow through the medium. This air flow was modeled by the 
continuity equation and Darcy's Law . The porous medium is 
characterized by a number of parameters including porosity , 
specific internal area (Z), tortuosity (7) and permeability 
(m). These parameters are defined as follows. Given a unit 


square cell of dimension D, and particle diameter d, we have 


pe 1 - « / 4 (d/D)- Giga) 
28 1/2 n d-/D- 6) 
Ee 1.4 (1.3) 
m - Sp?/z? (1.4) 


TSn p Mis defined as the ratio of void volume per 


m o икили The specifie internal area, Z, is the 


165 


ratio of internal surface area to pbull ps n ии. The 
tortuosity is defined as the ratio of the length of the 
flowpath of a fluid particle to the straight line distanci 

The particle's size decreases as the carbon is consumed 
and all geometric properties which depend on particle diam- 
eter are functions of time and position . This model 
considers that the carbon matrix remains rigid as the 
particle diameter decreases and thus porosity increases with 
combustion. 

The magnitude of REYNOLDS number defines the motion of 
the fluid which can be molecular , viscous or inertlal MN 
the case of porous media , the flow is dominated by viscous 
and inertial effects. Darcy's Law for the range of REYNOLD's 
number where viscous effects dominate is the equation that 
governs the fluid flow in porous media. Two velocities which 
distinguish the flow through the medium are filter velocity 


Q, and pore velocity u. Filter velocity is given by 
Q m mia EOXI) (1.5 


where u is the viscosity of the fluid and P is the pressure. 


The pore velocity u is given by 


U = Q/p (1506 


The hypothesis of the Dupuit-Forcheimer assumption is 
that the local pore velocity is greater than the filter 
velocity. The actual velocity in a pore is a functions 
position within the pore . The Dupuit-Forchermer rela tmo 
defines an average velocity in a pore. 

The model of combustion adopted was the model of N. N. 
Semenov (2). The relation of reaction rate to temperature 
and oxygen concentration, and the interaction of heat gener- 
ation and heat transfer are fundamental to this model . If 


the heat transfer dominates, extinction will oeecur 00057 
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heat generation dominates then combustion vill result. The 
reaction rate equation vas adopted from Arrhenius s Lav for 
a simple reaction vhere the rate depends on the concentra- 
tion of the reactants and not on the products. The heat 
generated by an exothermal reaction is obtained by multi- 
plying the reaction rate by the heat of combustion 

A study of the curve of heat generation versus tempera- 
ture presented by Frank-Kamenetskii distinguishes two phases 
during combustion. During the initial kinetic phase , the 
rate of reaction and the temperature are lowest. During This 
phase , there is an excess supply of oxygen. In the kinetic 
regime the  reation rate increases exponentially with 
increasing temperature. The kinetic regime is followed by 
the diffusion regime. This phase is distinguished by high 
temperatures and reaction rates , and the reaction is 
limited by a lack of oxygen. 

It is assumed that the chemical reaction for the combus- 
tion process produces carbon monoxide and carbon dioxide. 
The ratio of the mass rates of carbon monoxide to carbon 
dioxide depends on the temperature. Increasing temperature 
ик 4165 in an increase of this CO to CO» ratio. 

The model of heat transfer includes three basic mecha- 
nisms : conveeryon . conduction and radiation . The heat 
transfer equation used by VATIKIOTIS (1) also includes a 
heat generation term due to combustion. All thermophysical 
properties were treated as temperatures dependent. These 
properties include viscosity, conductivity and density. 
Energy balance equations were constructed for control 
volumes of the porous medium and the air flowing through the 
porous medium 

A third field equation of the model was obtained by a 
mass balance of oxygen . This equation includes molecular 
ENEIuSson. convective transport and a sink term to account 


for the depletion of oxygen due to combustion. 
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These equations resulted in a system of coupled, nonli- 
near transient field equation . These field equations, 
together vith boundary and initial conditions , define the 
problem. They vere solved numerically by using the Galerkin 
finite element method. Details of the formulation are 


presented in reference l 
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TU EEEGİ” OB HEAT INPUT ON COMBUSTION 


A. INTRODUCTION AND DESCRIPTION 

In this section the effect of heat input on the combus- 
tion behavior of porous medium was investigated. The porous 
medium is subjected to a heat source on its surface . Two 
parameters of main importance are considered: The magnitude 
of the heat flux of the source to the porous medium (SQ), 
and the time that the heat flux remains active (TQ). 

To do this analysis, all parameters of the system are 
fixed except the time period of heating (TQ) and the amount 
of heat flux (SQ). 

The purpose of this section is to determine the relation 
Between SQ and TO which will result in combustion. This is 
achieved as follows: 

For a specific value of heat flux SQ, the program is run 
ви Һат ЕГегепк values of TQ until combustion occurs at (TQ), 
2 extinction occurs at (10). , vhereö is a small time 


period . Here ve tookğ seconds 
m ТОдд 5 (CZ) 


Eur cases of heat flux were studied in this section: 


CASE II.1 SQ - 18000 BTU/ft.sq hr 
CASE II.2 SQ = 20000 BTU/ft.sq hr 
CASE aL. 3 5100-0000 BTU/it.sa hr 
CASE II.4 SQ = 40000 BTU/ft.sq hr 


The fixed parameters used in all cases are: 


Ambient temperature = 80 F 
Ambient temperature - 2117 lb/ft.sq 
Tortuosty = 1.400 


Filament diameter = 0.0004167 ft 
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Thickness of matrix laminate - 0.0004167 ft 
Thickness of porous medium - 0.02083 ft 

Gas constant for air = 5272271671: 1 
Conductivity of filament = 86 BTU/1bmiH sns 
Specific heat of filament = (О 70 5 Rlpm Ci 
Emissivity of filament - 0.90 

Shape factor for int.HF.XFER coefficient = 1.00 


Characteristc length of medium = 1 inch 
Heat of reaction = 14090 BYU/CE 
Reaction order = 0.50 

Stochiometric ratio (ftuel/alr) “7057 


Reaction coefficient - 2065000 Ibm, Gian 
Activation energy coefficient = 28840 deg R 


Pressure differential across thickness - 50 lb/ft.sq 
Initial carbon and air temperature = 80 F 
Initial uniform oxygen concentration = 0.0172 lbm/ PENES 


The equations used in this program permit three sets of 
boundary conditions that are approximations to phys m 
situations. The boundary conditions selected are the best 
approximations of a one dimensional model. The set of 
boundary conditions used in this case are typical of thermal 
flov reactor vithout radiation from the boundary surfaces 
for the problems of this investigation, the Eollowimsa 


boundary conditions vere used 


- No heat transfer from either end of the porous medium. 


di dx = 0) at "X/L.s.0 anda ии ди (2:20 
- Danckverts" boundary conditions : Conditions for the 
air temperature and oxygen concentration vere assumed. This 
means 
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Ka öTA / öx = PaCau(Ta- T.) XE 05) 
7:00 X/L 5 
poo /6x =u Coe $ + oo) X/L (2797 
R. —. 
FK.) 916 9£ = co E nö X/L (295 
Hed | ox = p-,€,U(T., - Te) DE AT, 00 d 
E gx = 0 XIL 6 
D.99/öx = ии фә) XIL 5 
57 PROCEDURE 
For each case , the amount of heat flux (SQ) vas fixed 


and several values of shut off time (TQ) vere chosen until 


2: of TQ. and TQ. were found such that 


io. 


n 


where 


2 


8 


nən 


IQ, * 8 


DS marl "mumber. 


EA 


710) 


(2501) 


Cee Ee 

The results for each case were obtained in both numer- 
ical and graphical form . The behavior of carbon tempera- 
ture, oxygen concentration and reaction rate for combustion 
and extinction are shown as surfaces versus time and posi- 
tion The program also calculates the air temperature in the 
porous medium. However the air temperature surfaces are not 
shown as they are quite similar , in most cases to carbon 
BSmpera unes: 

There are three surface graphics for each heat flux 
case: 

-Carbon temperature versus position and time 

- Oxygen concentration versus position and time 


-Reaction rate versus position and time 


In adition , there are two curves for each heat flux 
(SQ) , one showing the variation of carbon temperature 
versus time for combustion and the other showing the varia- 
tion of carbon temperature for extinct omik These is two 
curves are for the point x/L=0. (i.e. where the air enters 
the porous medium) 

l. “CASE 11-L: 50 13000; DT Ee so 

The values for shut off time bounding extinction and 


combustion are 


for combustion TQ, 46 seconds 


for extinction INO? 47 seconds 


During the transient analyses , system parameters 
varied . for each analysis , the maximum and minimum values 
of some of these parameters provide some insight into the 
character of the particular case . These parameters and 
their values are given in Table I 


The graphical results are shown in Figures 2.1 to 
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Figure 2.1 Temperature vs X/L and time 
Ber SO 76000 BTU/İt.sq hr and TQ = 46 seconds. 


CRRBON TOPE ERATURE if Pe 





Figure 2.2 Temperature vs X/L and time 
for SQ = 18000 BTU/ft.sq hr and TQ = 45 seconds. 


OXYGEN CONC. SURFACE FROM GREP SE 


OxrGEN CONCENTRATION (xJ10Q0CO! 





Figure 2 Oxygen concentration vs X/L and time 


for SO = 18000 BTU/ft-sqg nr andın. 46 seconds. 
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Figure 2.4 Oxygen concentration vs X/L and time 


for SQ = 18000 BTU/ft.sq hr and TO = 45 seconasa 
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Figure 2.5 Reaccion rate vs ^/L anc time 


for SQ =18000 BTU/ft.sq hr and TQ = 46 seconds. 
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Figure 2.6 Reaction rate vs X/L and time 


2717510 000 BIU ft. sq hr and TQ = 45 seconds. 
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2. CASE 11-2 SÜ — 20000 BTU FEC ip 
For this case the heat flux duration times for the 


conditions of extinction and combustYyoDW e 


35 seconds 


for combustion po 


for extinction TOF 34 seconds 


During the transient history of this case several 
parameters varied. The maximum and minimum of these changing 
variables are in the Table II . It is necessary to emphasize 
that these extrema results do not occur at the same time or 
same position (X/L). They are maximum and minimum values, 
obtained at different times and position during the tran- 
sient process. These values shov the characteristc behavior 
of this case 

The behavior of carbon temperature , oxygen concen- 
tration and reaction rate are shown in Figures 2.0 to) Za 
These surfaces are an important characteristic of this case 
because they show the evolution of the combustion process or 
extinction process for each position during the transient 
time 

2. (ASE II-5 50 = 39000 pTU/ fU sq hu 
The values of heating time that result in combustion 


and extinction for this value of heat f iita: 


for combustion lü. — 1535555: 
for extinction TQ, = 14 sseeonte 
During the transient analysis , some variation of 


values occured to several parameters. The limits of varia- 
tion , maximum and minimum values , are shown in Table III 
The process of combustion and extinction have Particulas 


characteristics that determine the variations i BE p Ди ии 
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Pieurew ao Temperature vs X/L and time 
for SQ = 20000 BIU/£t.sq Br and 107 0076 a 
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Figurer mi9 Temperature vs X/L and time 
for SQ - 20000 BTU/ft.sq hr and IQ € 3476660777 
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Psesure 2.10 Oxygen concentration vs X/L and time 
27::707- 0000 BTU/ ft. sq hr and TQ - 35 seconds. 
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Figure 2.11 Oxygen concentration vs X/L and time 
1057707 20000MBTU q hr and TQ = 34 seconds. 
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Page eee a Reaction rate vs X/L and time 


for SQ = 20000 BTU/ft.sq hr and TQ = 3556660907 
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Eigure 2 l3 Reaction rate vs X/L and time 


for SQ -20000 BTU/ft.sq hr and 10::02c6 60 
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FOR NUMBER POINT 1 AND SQ=20000 


Figure 2.14 Extinction and combustion carbon temperature 


for position be 0.. 
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of these parameters. The maximum and minimum values give a 
general idea of system behavior . These limiting values do 
not occur at same position (X/L) or same time 

A better idea of the development of temperature , 
oxygen concentration and reaction rate can be obtained from 
Figures 2.15 to 2.20 These surfaces shov the values of these 
parameters for each point (X/L) during the transient time 

Figure 2.14 shows that at X/L - 0. the transient 
behavior for extinction and combustion are the same until 
r. 14 seconds , the time the heat flux is terminated 
After this time the temperature for the extinction case 
decreases to ambient temperature ; while for the combustion 
case , the temperature increases reaching a peak temperature 
of 3200 F at 28 seconds and thereafter decreases to an equi- 
librium combustion temperature of about 2500 F after about 


50 seconds (see Figure 2.14). 


4. CASE II-4 SQ = 40000 BTU/ft.sq hr 
The heat flux of this case produced combustion and 


extinction of the porous medium at the following times 


for combustion IQ. - 9usseonss 
for extinct tom TQ. = 6 Second: 
During the transient process , the output votum 


program showed parameters th at didn’t change with time and 
others that did . The variation of these parameter's values 
show the development of the combustion and extinction 
phenomenon for this value of SQ . The maximum and minimum 
values of several parameters are shown in Table IV . Again 
these values do not correspond to the same time 

The characteristics of this case , carbon tempe: 


ture , oxygen concentration and reaction rate aSr fone Eon 
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Етрига 15 Temperature vs X/L and time 
for SQ = 30000 BIU/ft.sq hr and TO = 15 6686557 
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Figure 2 ho Temperature vs X/L and time 
for SQ = 30000 BTU/ft.sq hr and 10:70 010 5 0 7 
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Fipwne 2.17 Oxigen concentration vs X/L and time 
15:50 -30000 BTU/Tt.sq hr and TQ = 15 seconds. 
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Higure Z18 Oxigen concentration vs X/L and time 
25777-— 0000 ETU/ ftU.sq hr,and TQ = 14 seconds. 
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Figure Zeer? Reaction rate vs X/L and time 
for SQ - 30000 BTU/İt.sq hr and 10 45 secon 7 
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Figure 2.20 Reaction rate vs X/L and time 
for SQ = 30000 BTU/fE.sq hr and 10:-410:“55 = 
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of time and position (X/L) are shown as graphic surfacesmam 

Figures 2.22 to 2.27 . These graphical results give ane 

about the behavior of each point of the tEransitemeeni stem: 
For X/L = Q , both the combustion and3 extinc:s22 


cases show similar behavior to this time equal to about 8 


seconds , when the heat flux is terminated . After this 
time, for the extinction case , the temperature decreases to 
ambient temperature . For the combustion case , the temper- 


ature at X/L = O increased to about 3300 F , and thereafter 
decreased to around 2500 F (temperature of combustion) . 
Both the extinction and combustion cases achieve new equi- 
librium temperatures , about 2500 F for the combustion cases 
and ambient temperature for the extinction cases , after 35 


seconds (see Figures 2.28 ). 


тили 
Here , some observations are made about the effects of 


heat flux input and duration of heat flux on the combustusm 


process. 
l. Power Relation 
For each case , a pair of values (50, (100004. 
obtained . A plot of these point is shown in Figures 2i 
and 2.30 . It is observed that the relation between SQ and 


TQ of these graphics yields the approximate power relation 


(equation 2.12) 


2. Relation of Temperature and Oxygen Concent. See M 


During the initial heating of the porous medium , 
the temperature increases while the oxygen 158 being 
defleted. After the heat flux 1S removed , either extinction 
or combustion will proceed. If extinction occurs 77002 
carbon temperature decreases to ambient temperature while 
the oxygen is being restored to ambient concentreition 

If combustion occurs „then the oxygen inside the 


medium goes to zero and remains there . The only oxygen 
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Figure 2.22 Temperature vs X/L and time 
for SQ = 40000 BTU/ft.sq hr and TQ = 9 seconds. 
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Figure 2.23 Temperature vs X/L and time 
for SQ - 40000 BTUÜ/İt.sq hr and TÜ - бикә ии 
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Figure 2.24 Oxygen concentration vs X/L and time 
for SQ = 40000 BTU/ft.sq hr and TQ = 9 seconds. 
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Higure 2:25 Oxigen concentration vs X/L and time 
for SQ - 40000 BTU/ft.sq hr and TQ = 8 seconds. 
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Figure 2020 Reaction rate vs X/L and time 
for SQ - 40000 BTU/ft.sq hr and 1@ = Seconds; 
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Figure 2.27 Reaction rate vs X brand ии 
for SQ = 40000 BTU/ft.sq hr an3 TOP 656075 5 
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TQ. = 3.4 107 ини (2.129 


present is at the entrance at X/L = 0 . The temperature 
continues to increase to'a maximum , and thereafter the 
temperature decreases . Eventually the entire medium would 


achieve a uniform combustion temperature 


3. Equilibrium Temperature for Combustion 
In all combustion cases , the behavior was similar 
During the initial combustion period , the temperature of 


the medium at X/L = O rose between 2000 F (CASE II-1) to 
3400F (CASE II-4) . At this time , on the other side of the 


medium at X/L = 1. , the temperature of the medium was 
approximately 1500 F . This during thus early trans Temi 
stage there is aAT of 1300 F (for CASE II-l) and 1900 F 
(for CASE II-4) across the medium . As time proceeds Mime 
temperature at X/L = O0 decreased , while at X/L = l. , the 
temperature increased . That is , with progressing time , 
the medium moves towards a uniform temperature . In the 
present case , this equilibrium temperature of combustion is 


around 2400 F to 2500 F for all cases 

The mechanism of equilibrium of the temperature 
during combustion is due to the heat transfer Dy" air 797 
the air flows through the porous medium it transports the 
heat from the hot front part of the medium to the cooler 
back part , until all points of the medium have the same 


temperature 
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LLL BBrRECIOF DHECKNESS#ON GOMBUSTION 


A. INTRODUCTION 

The porous medium has several geometric parameters that 
can have direct influence on the combustion process. The 
effect of the thickness parameter on system behavior is 
analysed in this section. All other parameters are fixed 
except the thickness of the medium. 

Ihe purpose of this section is to determine the relation 
of the thickness on the minimum temperature to start combus- 
Eon (i.e., the ignition temperature). 

For a specified value of thickness, the program is run 
with different values of initial temperature. This initial 
temperature is equal at all points of the medium. The lowest 
value of initial temperature that results in combustion is 
B5”  enition temperature Tc, Extinction occurs ас Тс-д , 
vhereğ is a small increment of temperature. The value of6 
used in this parametric analysis vas 10 F 

Seven cases vith different thicknesses vere studied in 
sees: Section: 

ә ЛЛ THICKNESS = 0,25 inches 
CASE III-2 “THICKNESS = 0,50 inches 
GXSEEIDIT.SNENITORNESS = 0,7/5 inches 
CASE III-4 THICKNESS = 1.00 inches 
CAS Li SPBTHIGENESS = 2.00 inches 
CASE 111-6 THICKNESS = 4.00 inches 
CASERTI / THICKNESS =-6.00 inches 

merely Case 15 Characterised by dts thickness. The fixed 
parameters used in all cases are: 

- Ambient temperature = 80 F 

- Ambient pressure = 2117 lb/ft. sq. 

z ortuosity = 1.400 

- Filament diameter - 0.0004167 ft. 
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- Thickness of matrix laminate - 0.0004167 ft. 

- Thickness of porous medium 0050205 MD 

- Gas constant for air = 53:34 15777£ lone 

- Conductivity of filament = 86 BTU/l1bf HF 

- Specific heat of filament = 0.703 lbm/cf. 

- EMUSSIVIty Of filament. = Ue 

- Shape factor for int.HF.XTFER cocificiTenlz??”” 
- Heat of reaction = 14090 Blur ca 

- Reaction order = 0.50 

- Stoichiometric ratio (fub arp > 

- Reaction rate coefficient = 2065000 lbm/CF hr 

- Activation energy coefficient = 28.840 deg R 

- Pressure differential across thickness = 50 1b/ft? 


=- Initial uniform oxygen concentration 702700 lbm/ ft? 


Be PROCEDURE 

For each case, the characteristic length of the medium 
(thickness) was fixed and values of the uniform initial 
temperature were chosen until the value of temperature for 
combustion (T,) and temperature for extinction (150442. 


approximately equal (T. = T.) or (l 11 7” 


C RESULTS 

For each case, the results were obtained in both numer- 
ical and graphical form. The behavior of carbon temperatus 
oxygen concentration and reaction rate for combustione 
extinction are shown as surfaces versus time and position. 

The numerical results from the program show the varia- 
tion of several parameter values during the transient anal- 
ysis. The air temperature inside the porous medium, the 
griphite temperature and the oxygen concentration are 
presented at each position (X/L) at intervals of time” 0005 
system properties are also given . These properties charac- 
terize the process of combustion and extinction p ин 
different thicknesses. The maximum and minimum do not occur 


at the same point or at the same time. 
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EE ENEE in the form of surfaces. 
These surfaces show the behavior of carbon temperature, 
oxygen concentration and reaction rate versus time and posi- 
mö (X/L) for the combustion and extinction processes. 

l. CASE III-1l Thickness = 0.25 inches 

The values of the initial uniform carbon temperature 


bounding the extinction and combustion processes are: 


1100 F 
for ex inction Te = 1090 F 


tor combustions IC 


The development with time of the behavior of the 
graphite temperature, oxygen concentration and reaction rate 
during the combustion and extinction processes are shown in 
Figures 3.1 to 3.6 The transient parametric analysis provide 
some insight into the characterers of this particular case. 
The maximum and minimum values of these parameters are given 
in Table V 

2. CASE III-2 Thickness = 0.50 inches 
The combustion and extinction process, — case, 


start vith the folloving uniform carbon temperatures: 


960 F 


mor extinetion le 9502F 


for combustion Tc 


For this case, the graphical results for graphite 
temperature, oxygen concentration and reaction rate are 
shown in Figures 3.7 to 3.12 . The characteristics of the 
changing parameters in this case, are shown by the maximum 


“ud minimum values in Tables VT 
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Figure 3 l Temperature vs X/L and time 
for thickness = 022055 


Initial carbon temperature - 1100 F. 
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A 22 Temperature vs X/L and time 
for thickness ia 


Initial carbon temperature = 1090 F. 
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Бириси |) Oxygen concentration vs X/L and time 
for thickness = 0.25 in 


Initial carbon temperature = 1100 F. 
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Figure 3.4 Oxygen concentration vs X/L and time 
for thickness = 0.25 in 


Initial carbon temperature = 1090 F. 
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pomo 5 Reaction rate vs X/L and time 
for thickness = 0R omin 


Initial carbon temperature = 1100 F. 
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Figure seo Reaction rate vs X/L and time 
for thickness “230094 


Initial carbon temperature = 1090 EF. 
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Figure 3.8 Temperature vs X/L and time 
for thickness - 095 M 


Initial carbon temperature - 950 F. 
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Figure 3.9 Oxygen concentration vs X/L and time 
tor thickness = 0.50 in 
Initial carbon temperature - 960 F. 
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Figure 3.10 , Oxygen concentration vs X/L and time 
forşehickness - 0:50"1n 
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Initial carbon temperature = 950 F. 
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Бири ии Reaction rate VS X/L and time 
kor thickness - 0.501 
Ee carbon temperature - RE 
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DE rc ə Reaction rate VS X/L and time 
for thickness = 0.50 in 
Initial carbon temperature = 950 F. 
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3. CASE III-3 Thickness - 0-759219 NES 
During this case the minimum uniform temperature 
that starts the combustion process and the maximum uniform 


temperature for extinction are: 


for combustion Tc “7065 


for extinction Te $59 


During the transient analysis, the system parameters 
variedad. Por kurs case , the maximum and minimum values of 
some of these parameters provide some insight into the char- 
acter of the case. These parameters and their extreme values 
are given in Table VII 

Figures 3.13 to 3.18 show the graphical results 128 
carbon temperature, oxygen concentration and reaction rate 


fOr CHIS tease. 
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EEN Temperature vs X/L and time 
for thickHess “ 0075806: = 


Initial carbon temperature = 900 F. 
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TEMPERATURE SURFACE FROM GRAF3E 
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Temperature vs X/L ana time 


Figure 3.14 


(75 in 


for thickness 


GER, 


Initial carbon temperature 
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Oxygen concentration vs X/L and time 
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Figure 3.1 
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0 75 in 


Initial carbon temperature 


for thickness 


00 F. 
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OXYGEN CONC. SURFACE FROM GRAF SE 


8 
s 
E 
Í 
5 
5 
B 





Figure 3.16 Oxygen concentration vs X/L and tıme 
for thickness = 0.73 in 
Initial carbon temperature = 890 F. 
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Ereure mi; Reaction rate vs X/L and time 
for thickness = OM Drain 


Initial carbon temperature = 900 F. 
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Frere. Ig Reaction rate vs X/L and time 
for thickness - 0.75 in 


Initial carbon temperature - 890 F. 


4. CASE III-4 Thickness - 1.00 inches 
The program was run with several values of initial 
uniform graphite temperature . The temperatures bounding 


Semoustion and extinction are: 


Zor combussronulc = 850 F 


for extinction Le = 840 F 


The variable parameters of this case attain maximum 
and minimum values as given in Table VIII 

EM ens  c4S€ ehe Şersphical results for carbon 
temperature, oxygen concentration and reaction rate are 


55:10 in Figures 3.19 to 3.24 
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Elcuze” 19 Temperature vs X/L and time 
lor eluekness = 1.00 in 


fuma carbon temperatüre = 850 F. 
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Pigure 3.20 Temperature vs X/L and time 
nurun ıl 00Pin 


Initial carbon temperature - 840 F. 
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Figure 3r Oxygen concentration vs X/L and time 


for thickness = 1 00min 


Initial carbon temperature 658055: 
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Figure 3.22 Oxygen concentration vs X/L and time 
for thickness = 1.00 in 


Initial carbon temperature = 840 F. 
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Figure 3.23 Reaction rate vs X/L and time 
tor thickness = 1.00 in 


Initial carbon temperature = 850 F. 
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Fligüurex3 23 Reaction rate vs X/L and time 
osi chnesse- 1200 in 


Initial carbon temperature - 840 F. 


67 


66-91 GU De Zo LI 40°82 30ƏTƏOT77Ə02 UOL12297107 


8916070 801170 {cso O = Glee 0 xəqunu spTouzZəşq 
989 GEST SIZ SC81 (14/33) £at»o[sA^ 210d 
2921- DS 68cl- 09€ - (nə 37/q1) 3auərpe:43 Əzrnssərd 
9907 "b 1907 ZUM o (bS:33/q1) əasnssəud 
41971 "1811 1191 1161 (or-01)(bs*33) Азтттдеәиләда 
unurTuTu: unurxeul unuTuTu — Unurxeul 
uorjsnquioo2 иотдоитдхә S193o9ue1ed 


S4HƏNİ 001 - SS3NXƏİIHL WO4 SS32ONd AHL INIANA SAALINVAVA G314VA 
IIIA 38V1 


68 


EE, SEET res 
The combustion and extinction process, in this case, 
are characterized by the following uniform initial graphite 


EE CO re SS 


for combustion Ic = 760 F 


for extern tree 750 F 


The graphite temperature and other parameters varied 
during these processes . The variation of these parameters 
characterize the system. The maximum and minimum values of 
these parameters are in Table IX 

The graphical results for carbon temperature, oxygen 


concentration and reaction rate are shown in Figures 3.25 to 
70 
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Figure 25 Temperature vs X/L and time 
"r ehiceness = 2:008in 


Initial carbon temperature - 760 F. 
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Figure 3.26 Temperature vs X/L and time 
for thickness - 2.00 in 


Initial carbon temperature = 750 F. 
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Figure 3.27 Oxygen concentration vs X/L and time 


for thickness - 2.00 in 


Initial carbon temperature = 760 F. 
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OXYGEN CONC. SURFRCE FROM GRRF3E 
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Figure 2:25 Oxygen concentration vs X/L and time 
for thickness - 2.00 in 
Initial carbon temperature - 750 F. 
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Fıgur2 3.29 Reaction rate vs X/L and time 
for thickness - 2.00 in 


Initial carbon temperature - 760 F. 
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Вирите "ир Reaction rate vs X/L and time 
for thickness - 2.00 in 


Initial carbon temperature = 750 E, 


6. CASE III-6 Thickness = 4.00 inches 
The minimum value of the initial graphite tempera- 
ture to start combustion of the porous medium and 2 
maximum value of the initial graphite temperature that 


results in extinctlon for this care aie.: 


for combustion Ice p E 


for extinction le = 605 


Other variable parameters of this case are shown in 
Table X . The values of these parameters are characteristic 


of the combustion and extinction processes for this value of 


thickness. 
In this case, the graphical results “huvə 


behavior of carbon temperature , oxygen concentration and 


reaction rate in Figures 3) 3 o DX 
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Temperature vs X/L and time 


Figure 90201 


2005 


Initial carbon temperature 2760-7017 


for thickness 





Temperature vs X/L and time 


Figure 3-32 


And 


Initial carbon temperature - 


for thickness 
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P3eurew 999 Oxygen concentration vs X Bande time 
for thickness - 4.00 in 


Initial carbon temperature - 690 F. 
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Figure 3.34 Oxygen concentration vs X/L and time 
for thickness = 4.00 in E 


Initial carbon temperature - 680 F. 
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HS la Reaction rate vs X/L and time 


for thickness = 7700030 


Initial.carbon temperature.- 690 F. 
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Figure 3.36 Reaction rate vs X/L and time 
for thickness, - 40D 


Initial carbon temperature = 9§6800F. 
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7. CASE III-7 Thickness “ 6:00:1m5155 
Below is given the initial uniform carbon tempera- 


ture for combustion and ex ri r ion: 


650 seconds 
640 seconds 


for combustron Ne 


for extinceon272 


The variable parameters for this case are given in 
Table XI 

The Figures 3.37 to 3.42 show the development, m 
carbon temperature, oxygen concentration and reaction rate 
at each position during the combustion and EEGEN 


processes. 





Figure 3.91 Temperature vs X/L and time 
for thickness = 6.00 in 


Initial carbon temperature = 650 F. 
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Figure 


Oxygen concentration vs X/L and time 


6.00 in 


Initial carbon temperature 


for thickness 


SUR. 
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Figure 3.40 Oxygen concentration vs X/L and time 


for talckness 85 00” 





Figure 3.41 Reaction rate vs X/L andıran 
for thickness = 6.00 in 


Initial carbon temperare cL 
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Figure 3.42 Reaction rate vs X/L and time 
for thickness - 6.00 in 


Initial carbon temperature - 640 F. 


De SUMMARY 

Here, some observations are made about the effects of 
thickness on the combustion process. 

l. Power Relation 

For each case a pair of values (thickness, combus- 
tion temperature) was obtained. Four plots of this data are 
shown in Figures 3.43 to 3.46 

Figures 3.43 and 3.44 show thickness versus combus- 
tion temperature , and Figures 3.45 and 3.46 show thickness 
versus the time it takes for the oxygen concentration at X/L 
- l. to reach. 

It is observed that the relation between thickness 
and combustion temperature , figure 3.43 has a form which 
approximates an exponentially decreasing function. 

The log-log plot ( Figure 3.44 ) is close to a 
linear form, which yields the approximate power relation. 

The temperature is in degrees fahrenheit and thick- 


mess is in inches 
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T = 850 L^0,16555 (AN 


A EOMPaArison Of the “Computer program S results and 
the formula's results are given in Table III. The maximum 
difference between the analytical results and the power 
equation results is less than 3 percent over the range of 
thickness between 0.25 and 6.00 inches 

This shows that the formula is a good approximation 
for these cases. 

2. Combustion Start Time 

During the initial phase of the combustion process, 
the oxygen concentration at all points of the porous medium 
continuously decreases to zero. for all cases, the point of 
the porous medium that first reached zero was X/L = 1. The 
time it taken for that the oxygen concentration at X/L = l. 
to reach zero is important because it is related to the time 
Bnat combustion starts. 

Figure 3.45 , shows that this time does not increase 
linearly with the thickness. The log-log plot of these 
values, Figure 3.46 , shows that this log-log plot is 
approximately a straight line. 


The equation for this line is: 
T* = 226.8 L1,683 (522) 


This equation can be used to give an estimate for 


ene time it takes for combustion to start. 
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THICKNESS VERSUS TEMPERATURE 
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Thickness versus combustion temperature 
Retangular plot. 
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Figure 3.44 Thickness versus combustion temperature 
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Figure 3.45 Thickness vs time for oxygen concentration 
at X/L = F to ee 


Retangular plot. 
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Figure 3.46 Thickness vs time for oxygen concentration 


reach zero log-log plot. 
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A. INTRODUCTION 

The permeability or hydraulic conductivity of a porous 
medium is the capacity of air flow through the porous 
medium. The air provides the oxygen for combustion . On the 
other hand , air flow provides for convective heat transfer 
within the medium . Thus air flow provides for heat genera- 
tion from combustion , as well as heat transfer by convec- 
tion . The interaction of these two mechanisms determines 
Enether combustion or extinction will occur . If the heat 
generation of combustion is greater than the heat transfer 
of convection , then combustion will occur . If the heat 
transfer mechanism dominates , then extinction will occur 


The equation used for permeability in the program is 


m - 0.2 p? / z” (4.1) 
Vhere m is the permeability , p is the porosity and z is the 
specific internal area . The expressions for porosity and 


specific internal area are 


m (4/0) (4.2) 
EE 102 — (4/02) (4.3) 
z - 1/2 nz (42102) 75 


Bue equation 4.3 is for spherical particles and equation 4.4 


2: “lindrical particles . Equation 4.1 is the Kozeny 
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equation for permeability of porous medium . This equate 


fails for highly porous fibrous medik 





Figure 4.1 Geometry of a typical cell. 


This equation fails for highly porous fibrous mea. 

The tortuosity ,z , 1S the ratio of the lengeh oia 
flow path for a fluid particle to the sStrargiıu ğı 
distance. The tortuosity depends on the ratio d/D Inm 
Program tortuosity was assumed to have a constant value of 
1.4 

Normally the interior geometry of a porous medium can be 
inregular and complex . In such cases the internal struc mE 
is difficult to discribe analytically . The program used 
this investigation models the porous medium as consisting of 
cylindrical fibers or spherical particles in a regular poms 
odic geometry as shown in Figures 4.1 

During the combustion process carbon is beine “onur. 
ously consumed and as a result, the internal geometry 1S 
changing. The medium becomes more porous as the carbon 
decreases. Therefore, all properties which depend upon fiber 


or particle diameter are functions of time and posicione 


a 
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model used in the program assumed that the carbon matrix 
remains rigid as the particles diameter decreases , and thus 


, the porosity increases with combustion 


ES C PROCEDURE 


To study the effects of permeability on combustion , all 


parameters of the program were maintained constant , except 
the permeability . In order to change the permeability and 
keep the porosity constant , it is necessary to keep the 


(d/D) ratio constant as seen in the equations 4.5 and 4.6 
This is accomplished by setting d/D equal tog . Then the 
numerator of equation 4.1 remains constant , however the 


denominator changes because d =q@ D and Z gives 


E eeneg d/D”) = constant / D (205) 
2 2 a2/n3Y- : 

Z = 1/2 gla d-/D~ )= constant / D (ano) 
Thus , the dimensions of d and D were changed propor- 


tionally so that the initial ratio (d/D) remained constant 
in all cases . Thus the porosity remained constant while the 
permeability was different for each case 


il ocoses cEUĞT€SC in this oseffion are 


CASE IV-1 PERMEABILITY = 0.00015 ft? (d-D-0.0025 in) 
CASE IV-2 PERMEABILITY = 0.00058 ft^ (d-D-0.0050 in) 
CASE IV-3 PERMEABILITY - 0.00233 ft? (d-D-0.0100 in) 
CASE IV-4 PERMEABILITY - 0.00933 ft? (d-D-0.0200 in) 


For each values of permeability the program was run with 
different uniform carbon temperatures for all points (X/L) 
This procedure was repeated until we obtaind the minimum 


temperature that results in combustion . This temperature 


om 


is called the combustion temperature (T_) | Extimomuem 


occurs at temperatures below T, 
0 (4.7) 


In this section this difference (08 ) was chosen equals 
10 F then 


Tr Oe) (4.8) 


Gas RESULTS 


For each case , the program provided the values of the 
carbon and air temperatures inside the porous medium , as 
well as other parameters which varied during the combustion 
process 

1. CASE IV-l Permeability = 0.00015 £t* (d=D=0.0025 in) 

For this value of permeability , the maximum carbon 
temperature that does not result in combustion (extinction 
temperature) and the minimum temperature that results in 


combustion (combustion temperature) are 


HH 
ll 


760 F 
7^5 0 NIS 


for combustion 


m. 
Y 


fOr Extinction 


The graphical results are shown in Figures 57777 
4.7 . These graphical surfaces show the carbon temperature , 
oxygen concentration and reaction rate at each position and 
time . The form of these surfaces shows the different 
behavior of these three parameters during the extinction and 


combustion processes. 
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Figure 4.2 Temperature vs X/L and time 
for permeability - 0.00015 ft? 


initial carbon temperatiure = 760 F. 
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Figure 4.3 Temperature vs X/L and time 
for permeability - 0.00015 ft? 


initial carbon temperatiure - 750 F. 
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Figure 4.4 Oxygen concentration vs X/L and time 
for permeability = 0.00015 ft? 
initial carbon temperatiure - /60 F. 
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Figure? 5 Oxygen comcentration vs X/L and time 
for permeability 009009 fr” s 
initial carbon temperatiıne7 8 
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Figure 4.6 Reaction rate vs X/L and time 
fer permeability = 0.00015 £t? 


initial carbon temperatiure - 7777. 
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Figure 4.7 Reaction rate vs X/L and time 
for permeability = 0.00015 ft” 


initial carbon temperatiure - 750 F. 
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2. CASE IV-2 Permeability = 0.00058 ft* (d=D=0.0050 in) 


Here the value of permeability is 0.00058 ft^ which 
results from a pore diameter of 0.0050 in . For this s 


the combustion and extinction temperatures are 


for combustion T. ~ EOM) 

for extinction DUM CE 
From the numerical output of the program , three 
graphical surface s were created ; the carbon temperature , 


oxygen concentration and reaction rate as functions or m 
and position . These surfaces provide information of the 
development of these three dependent variables during the 
combustion and extinction process . These surfaces are shovn 


on Figures 4.8 Co ATi 
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Figure 4.8 Temperature vs X/L and time 
for permeability = 0.00058 ft? 


initial carbon temperatiure - 850 F. 
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Figure 4.9 Temperature vs X/L and time 
7777557 0 00055 £€” 


initial carbon temperatiure = 840 F. 
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Figure 4.10 Oxygen concentration vs aa E me 
Бе J U 0 00058. FrZ 


initial carbon temperatiure = 850 F. 
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Figure 4.11 Oxygen comcentration vs X/L and time 
for permeability - 0.00058 ft? 


initial carbon temperatiure - 840 F. 
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Figure 4.12 Reaction rate vs X/L and time 
for permeability = 0.00058 ft? 


initial carbon temperatiure - 850 F. 
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Figure 4.13 Reaction rate vs X/L and time 
for permeability - 0.00058 ft” 


initial carbon temperatiure - 840 F. 


3. CASE IV-3 Permeability = 0.00233 ft? (d-D-0.0100 in) 
For this case the initial permeability is 0.00233 
ft? . With this value of permeability the combustion and 


extinction temperatures are 


0800 FE 
770071 


for combustion T 


Tor exei nction E 


The variation of the dependent variables during the 
combustion and extinction cases is shown by a numerical 
output of the program . These results are also shown in 
three graphical surfaces ; carbon temperature , oxygen 
concentration and reaction rate on each position (X/L) and 


time . These graphical results are shovn in Figures 4.14 to 
4.19 
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Figure 4.14 Temperature vs X/L and time 
for permeability = 0.00233 ft? 
initial carbon temperatiure = 980 F. 





Figure 4.15 Temperature vs X/L and time 
for permeability - 0.00233 ft? 
initial carbon temperatiure 7777070807 
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Figure 4.16 Oxygen concentration vs X/L and time 
for permeability - 0.00233 ft” 


initial carbon temperatiure = 980 F. 
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Figure 4.17 Oxygen comcentration vs X/L and time 
for permeability - 0.00233 ft” 
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Figure 4.18 Reaction rate vs X/L and time 
for permeability = 0.00233 ft” 


initial carbon temperatiure - OS UNES 
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Figure 4.19 Reaction rate vs X/L and time 
for permeability == OMG | ft” 


initial carbon temperatiure “ 970E 


To 


A. CASE IV-4 Permeability - 0.00933 ft? (d-D-0.0200 in) 


The combustion and extinction temperatures for this 
2 


Bec with initial permeability equal to 0.00933 ft” , are: 
tor comDustion [EDDIE 
ftor.exbinetron 27::155 1 

for this case , it vas not possible to obtain the graphical 


surfaces due to limitations of the graphical program 


D. SUMMARY 

Here , some observations about the effects of perme- 
ability on combustion can be understood by looking at all 
the results at one time . Table XIII contains the results of 
all cases 


The Table lists the temperatures for extinction and 


combustion ,the permeability (m) , the dimension of the unit 
cell (d=D) and time t* . The t* time is the time at which 
the oxygen concentration at a point reaches zerosvtirst 


That point of the porous medium that first reaches zero 
oxygen concentration in all these cases is the the position 
point X/L = 1 
l. Power Relation 

The pair of values (permeability , combustion 
temperature) for each case was plotted on a cartesian graph. 
5207 observation of this curve shows that the combustion 
temperature increases with increasing permeability. Figure 
4.20 shows that the combustion temperature increases mono- 
tonically and approaches an assindoticptate as permeability 
increases . The plot of these same points on a log-log 
570207  “Figures 4.21 , shows that its form is not quite a 
Straight line 

It is observed that the relation between combustion 
BEuperature (I-J) and permeability (m) yields the approximate 


power relation (equation 4.9) 
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75:15:57: (4.9) 

A comparation of the values of the combustion 
temperature obtained by the program and by the pover rela- 
tion is presented in Table XIV . This Table shovs that the 
difference betveen these results fail betveen 1.2 to 3.2 
percent . This result shovs that the above pover relation is 
a good approximation to estimate the temperature vhich vill 
start combustion for a given permeability 

2. Combustion Speed | 

In this section , the question of how permeability 
affects the time it takes for combustion to begin is 
addressed . We meassure the rapidilly of the combustion 
process by the t* time previously defined . Recall that t* 
is the time at which the oxygen first reaches zero at some 
point in the porous medium . In all cases the position X/L - 


l is the first point to achieve zero oxygen concentration 


The plots of permeability (m) versus t* , Figures 4.22 and 
4.23 , show that t* decays exponentially with increasing 
permeability and approaches an assymptotic value of t” = О 


as m goes to infinity 

À very interesting observation is obtained from the 
m ules of Table XIII . It is noted that the m;,4/m; (for i 
BBB Ü mains pratically constant (3.84 , 4.0L , 4.00) 
75: :16 U”)/ÜE”) ratio takes on the variable values 
1 01) :1/1391 , 1/2.65 , 1/3.56) 


An attempt vas made to obtain an additional run vith 
the permeability equal to 0.03733 ft? , hovever the combus- 
tion activity vas too rapid to achieve meaningful results 
Combustion vas achieved so rapidly that the gradient vere 
too large for the program to handle 

The permeability is a parameter that has great 


influence on the combustion process for porous medium 
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Figure 4.20 Permeability versus combustion temperature 


retangular plot. 
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Figure 4.21 Permeability versus combustion temperature 


log log plot. 
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Figure 4.22 Permeability versus t* 


Retangular plot. 
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Figure 4.23 Permeability versus t* 
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V. REACTION ORDER 


A. INTRODUCTION 
The expression for reaction rate for carbon reacting in 
air used in the program vas the expression resulting İrom 


the experimental data of Parker and Hotel 


n 








Hz A o, ex - 3 
c m ə = (5.1) 
C c 


Reaction order is the exponent of the partial pressure 
of oxygen in the expression of the reaction rate . The 
program permits the user to select any value for this param- 
eter . Im all of the previous chapters , a value of n equal 
to 0.5 was used . In this chapter we investigate the effect 


of n on the combustion process 


Eee) PROCEDURE 

In this section , all parameters are fixed to their 
basic values and only the reaction order is changed in each 
case 

The program was run , for each value of reaction order , 
with different value of the uniform initial carbon tempera- 
ture until we obtained the minimum temperature which results 
in combustion . This temperature is the combustion tempera- 
2:5: 1.) . The maximum temperature that results in extinc- 
2 is the extinetion temperature (T.). they differ by an 


dulinitesimalö (1, = To 0) . The magnitude ofó was chosen 


€ 
equal to 10 F for this study 

In order to study the influence of the reaction order on 
the combustion process , four different values of this 


parameter were chosen . These cases are 


ını 


CASEIN I Reaction Order 308975 


EE Reaction Order (1) 
CASE V-3 Reaction Werder” Фи: 
CASE V-4 Reaction Order = Vue 


The fixed parameters used in all cases are 


Ambient temperature = 80 F 
Ambient pressure - 2117 lb/ft.sq 
Tortuosty = 1.400 

Filamento diameter = 0.0004167 ft 


Thickness of matrix laminette = 0.0004167 ft 
Thickness of porous medium = 0.02083 ft 

Gas constant for ait =) 55404 IDF Aa 
Conductivity of filament = $6 BU Abner. уи 
Specific heat of filament = 0.703 lbm/Cf 
Emissivity of filament = 0.90 


Shape factor for internal HFxFER =1.00 


Characteristc length of medium = 1 inchess 
Heat of reaction. = “14090 EIER 
Stochliometric ratio (fuel/arr) = 0 25 


Reaction coefficient “””706252000 lbm, 6-77 
Activetion energy coefficient - 28840 deg R 
Pressure differential across thickness - 50 lb/ft.sq 


Initial uniform oxygen concentration 420 MM 1bm/ ft? 


C S RESULTS 

For each cases the results were obtained in both numer- 
ical and graphical form . These results show the ҹагтарины 
of some parameters during the combustion and extinction 
processes . The results obtained for each value of reaction 
order are presented below 

1. | CASE V-l Reaction Order TnT 

With this value of reaction order , the combust 

and extinction processes define the following characteristic 


Lempenatune- 


Lü 


980 o.F 
970 o.F 


POr COMBUSTION Te 


HƏ ne Edo Ta 


The analysis of the numerical output of the program 
shows that some parameters varied during the combustion and 
extinction processes . The limits of these parameters show 
BH ce character Jof the "Particular kind of process which 
occured 

The graphical results show three surfaces that 
represent the evolution of the carbon temperature , oxygen 
2: sretsion and rescelon ratesat each position during the 
transient history . These graphical surfaces are shown in 


55cures 5.1 to 5.6 
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Figure 5.1 Temperature vs X/L and time 
For reaction order = 0.25 
Initial carbon temperature = 980 F. 
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TEMPERATURE SURFACE FROM GRRF3E 





Figure 5.2 Temperature vs X/L and time 
for reaction order = 0.5 
Initial carbon temperature = 970 F. 
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Figure o3 Oxygen concentration vs X/L and time 
for reaction order =mi S ü 
Initial carbon temperature - 980 F. 
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Figure 5.4 Oxygen concentration vs X/L and time 
for reaction order - 0.25 
Initial carbon temperature = 9/0 F. 
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Figure 5.5 Reaction rate vs X/L and time 
for reactiem order = 0.25 i 
Initial carbon temperature - 980 F. 
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ЕЈјвиле 6 Reaction rate vs X/L and time 
for reaction order = 0.25 


Initial carbon temperature = 9/70 F. 


2. CASE V-2 2 Reaction, iu 
For reaction order equal to 0.50 , the inic mmi 
uniform carbon temperatures that separate the combustion and 


extınction processes are 


980 F 
for extinction 13 = 070 


for combust sop ra 


The program s numerical output included the values 
of some parameters that varied during the transient process. 
The variation of the values characterize the combustion or 
extinction process. The graphical rəsults are shovn in 
three surfaces: carbon temperature, oxygen concentration and 
reaction rate versus position (X/L) and time These graph- 


_— 


ical surfaces are shown in figures > mee m 
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Fipure 5.7 Temperature vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 850 F. 
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Figure 5.8 Temperature vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 840 F. 
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A a. 7 Oxygen concentration vs i/L and time 
for reaction order — 0m” 
Initial carbon temperature - 850 F. 
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Figure 5.10 Oxygen concentration vs X/L and time 
for reaction Order “05-05 
Initial carbon temperature - 840 F. 
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Pigūre 5.11 Reaction rate vs X/L and time 
for reaction order = 0.50 
Initial carbon temperature = 850 F. 
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Fipune 5.12 Reaction rate vs X/L and time 
for reaction order - 0.50 
Initial carbon temperature - 840 F. 
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3. CASE V-3 Reaction Order - gon 
The reaction order was set equal to 0./5 for this 


case. The combustion and extinction temperatures are: 


750 F 
740 F 


for combüustio mna 


for extincion Me 


A visualization of the development of the carbon 
temperature, oxygen concentration and reaction rate at each 
position and time are shown in the graphical surfaces of 


FPFigures@5 10:662271 
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Elcurell Temperature vs X/L and time 
for reaction order = 0.75 


Initial carbon temperature = 750 F. 
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Figure 5.14 Temperature vs X/L and time 
tor reac tion order = 0,75 


Initial carbon temperature = 740 F. 
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lrpure 5.15 Oxygen concentration vs X/L and time 
for reaction order = 0.75 
Ee EE SE - 750 F. 
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Figure 90016 Oxygen concentration vs Reeg 


for reaction order - 0.75 
Initial carbon temperature - 740 F. 
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Empeure E) Reaction rate vs X/L and time 
for reaction order = 0.75 
Initial carbon temperature = 750 F. 
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Eppure 5-15 Reaction rate vs X/L and time 
for *eaction order = 0.75 
Initial carbon temperature = 740 F. 


777057: Reaction Order = 1.00 


For this case, the temperature of the combustion and 


extinction processes are: 


for combustion Ta -u070. E 


for -extinction e: = 16600" E 


The Figures 5.19 to 5.24 show the behavior of the 
carbon temperature, oxygen concentration and reaction at 


each position and time. 


b. SUMMARY 
1. Power Relation 
For each case, a pair of values (reaction order, 


combustion temperature) was obtained.The points were plotted 


Me artesian coordinates in Figure 5.25 . A log-log plot is 
shown in Figure 5.26 , and a semi log plot is shown in 
Figure 5.27 . It is observed that the relation between 
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Figure 5.19 Temperature vs X/L and time 
for reaction order - 1.00 


Initial carbon temperature - 670 F. 


reaction order (n) and combustion temperature (T -) of topo 
5.27 yields the approximate power relation (equation 777 
Table 5.25 shows a resume of each case, combustion 
temperature, extinction temperature and the value of the 
combustion temperature obtained from Equation 5.2 and the 
percentage error . An observation of these results leads to 
the conclusion that this power relation is a good estimation 


for combustion temperature. 
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eege 5.20 Temperature vs X/L and time 
for reaction order = 1.00 


Initial carbon temperature = 660 F. 
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Figure 5.21 Oxygen concentration vs X/L and time 
for reaction order = 1.00 i 


Initial carbon temperature - 670 F. 
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Figure 5.22 Oxygen concentration vs X/L and time 
for reaction order = 1.00 
Initial carbon temperature = 660 F. 
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Figure 5.23 Reaction rate vs X/L and time 
for reaction rate= 1.00 


Initial carbon temperature = 670 F. 
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Figure 5.24 Reaction rate vs X/L and time 
r resctionlorder - 1.00 


Initial carbon temperature - 660 F. 
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TABLE XV 


COMBUSTION AND EXTINCTION TEMPERATURE 
FOR REACTION ORDER CASES 


temperate o 


extinction combustion A erron 
reaction program program power 
order output output relation 
2 770 980 0 ; 
50 840 850 863 Le 
(115 740 750 760 1.4 
E00 670 670 670 0.0 
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Peewee 5,25 Reaction order versus temperature 


Retangular plot. 
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REACTION ORDER VERSUS TEMPERATURE 
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Figure 5.26 Reaction order versus temperature 


Log log plor: 
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REACTION ORDER VERSUS TEMPERATURE 


COMBUSTION TEMPERATURE (F) 
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Figure 5.27 
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Reaction order versus temperature 


Semi log plot. 
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VE: CONSLUSTONS 


In each of the previous chapters, the effects of a 
single parameter on the combustion behavior of a porous 
medium was investigated. In Chapter 2, we considered the 
effect of heat flux (SQ) and its duration (TONOS 
combustion. In Chapter 3, the effect of medium thickness (L) 
on behavior was investigated. Chapter 4 was a study of the 
effect of permeability (m) on system behavior; and finally 
in Chapter 5, the effect of reaction order (n) was 
investigated. 

The major conclusions from these investigations are: 

- The relation between SQ and TQ is not linear , a power 

relation gives rasonable approximation. 

- During the combustion process, for the SQ cases, 
initially there is a large temperature gradient across 
the medium with T ( X/L = 0.) >> T ( X/L-1./77777 79 
passing time, this gradient decreases until there is a 
uniform combustion temperature over the entire medium. 
When this occurs, there 1S no oxygen within the 
medium, and the medium is continuosly consumed by 
Conbustiomn. 

- The relation between thickness and combustion tempera 
ture is not linear . Increasing the thickness results 
in decreasing of combustion temperature. If the 
thickness increased with a constant rate c (COS 
c ) the time that marks the start of combusti ici K 


constant , it increases with increasing thfotnes:l 


DA 
^ 


ires. t SPINA - b(L) ) where b depends upon the 
thickness 

- The relation between permeability and minimum combus 
tion temperature is not linear. 

- An increase of the permeability shows that the combus 


tion temperature increases to an assyinpliobrc sn 


1. 


- Permeability has a significant influence on the speed 
of the combustion process. When the permeability is 
increased, the time that marks the start of the 
combustion process decreases. For large values of 
permeability, the : combustion may occur in a fraction 
of a second. 

- If the permeability is increased with a constant rate 

TIME Di n time that measures the start ol 

combustion decreases 


c( m 


- The relation between reaction order (n) and combustion 
temperature (Tc) is not linear. 

- Àn observation of the graphical surfaces, for all 
cases of this investigation , shows a similar form for 
all combustion and extinction cases. borea 
combustion cases, it was noted that the oxygen within 
the porous medium being consumed with passing time. 
The first point at which the oxygen is depleted is 
always X/L = 1. With passing time the point with 
Sia reach zero oxygen, until all points have 
zero oxygen. At this time the only oxygen is at the 
interface of the porous medium and the ambient 
environment. At this time the combustion process 
becomes a surface recession phenomena. At the same 
time the oxygen is being depleted the temperature is 
increasing. Eventually when all the interior oxygen 
has been consumed, there is a temperature gradient 
aeross the medium with T.(X/L-0) > T (X/L=1). 

This thesis concludes with two recomendations . First 


5 


it is recommended that additional parameter studies he 


conducted . For example , it would be useful to determine 
the effects of pressure differential , conductivity , and 
other thermophysical properties on combustion . Another 


recommendation is to attempt to obtain a relation between 


combustion temperature and several parameters 


195 


LIST OF see ene eis 


Vatikiotis , C. S., A Combustion and Heat Tranfer Model 


for Porous Media , Phd.Dissertation , Naval 

Posgraduate School , Monterey , California , 1982 
semenov , N. N. , Chemical hinetics “and Cha Hec a 
Clarendon Press , Oxford 50” 


134 


mili 5505 EKLEULLON LIST 


No. Copies 


Departamento de Pesquisa e Desenvolvimento 2 
Ministerio da Aeronautica 
Brasilia , Brazil 


Centro Tecnico Aeroespacial 2 
Ministerio da Aeronautica f 
sao Jose dos Campos , SP , Brazil 12200 


Manor Antonio E. S. Serapiao , Brazilien air force 3 
Centro Tecnico Aeroespacial 
900106 710 «Campos , SF , Brazil 12200 


Associate Professor David Salinas , Code 69Zc 3 
Department of Mechanical Engineering 

Naval Posgraduate School 

Monterey ; California 93940 


porary , Code 0142 1 
Naval Posgraduate School 
Monterey , California 93943 


Departament Chairman , Code 69 | 1 
Departament of Mechanıcal Engineering 

Naval Posgraduate School 

Monterey , California 93943 


Defense Technical Information Center 2 
Cameron Station ` 
Alexandria , Virginia 22304-6145 


155 


















E da B rfr we ww mq vee q w [ 
| "www AA a e NN e A e 
- “nür 











M məv vv ” — 
a ? "ә: ДЕ. 
Ы = END o al 
7 Fs Tas qui Cd dad o 2 7 Кан A A A TA A Da e 
> a. 5 H ә k ve . 
9 "y 66m n 2111 rper ¡ba ne% La: y r.) KN KC p Ar P4 BBB as a e TTT TT 
"ə. 224 с CA Дад. 
| Ad da e ә ыы E A 
Š DH í | : 
1 e 
1 = 
= ° | 
a 
= > б 
E? LÀ en. . 
E . Ë j 
ES beier : 
0. pa R 
PX lët, J 
g : i 1 ¿ > - — = a 
- D - T 1 | - 
E E - a ҹ | d 
LJ - - = p Ы 
_ .. Ы - k — L. We == f : 
ber LE KE d). ies à : 
- E nə ~ wt Sen a E | 
P ae . -- n | | i 
LJ Q a basi ~n “ 1 we > z ! 
s b e ız = [€x E = “us hl - - Pod D 
fr um ES = y | 2 = 
= 2 ee — | “s i e " es = və ya ^ s Le PN - 2. 
- ES ə. ç “r. 2 F I j | ; : 
— ean də. vm E p = x" | : , ! V : 
LT > - aaa "o WK wv H = S - E PM te weng 2 ` 
: AM d E A 2.» o T Ы < ge nd d H a e: - Lodi s. ni adus 
4 T = — Pd "T Y “ər x“ WW ə = : E 
= amar AE EI 1... ECK Ae eg KR Ç : = 2 = 
-< ҝӱҹ Sed EEN E o e ә” aA s ` "LI. o ә/ ae ə 
= R => s - pr ду — M zn a ~ em is - = KH i. A nid - A e rore = ар —— o a q TES — SS ” d P sad 
op, a , . "E TNI. ~ e D LEX nə “dü... E - "T — m p. ms. if adi ndi sees ҹ̧ә ә 
D e. a La meng, dn. zg Zeg "hiem al » meie] g z a y "Y : : d | | zi I 
— ~ ә Рә q^ KSE CH ER E ~ “ә an - E a .. a de ә ә PLI keete c a eg y 
LN = D E Lu - ә ә бә e : = | : = n 
Wë do ip A PU a e ə. Аҹә I es ° - fms EI D n * = e ah a» d^ ^ we P - E T UU T. 
~ E . z N SE AE A ao E PT ~ A —p -. * PITE ... ZE [Jd 
EN - - m w - `. > | - e P a —ÀÓ > * v» fo c.m mt P um. Á— at kW 
2 ти - +. el XL e eng p - e E E E e S Lj - d - P go tv dani auam - ep 
m 2 < wa a EH past F x mem > - - e up tege" = Y 
= ҹа as m > 5.1... w ` v - A» a e EET - - - ya ~ P ww 20 6 ad — - ə - Рә. 
ə pm рә лт D = kb PIT H PT idi UN d .. E ke Kg = Ze - P ei "em o CE r ee e po 
= deeg wir erh deefe ect e ESA = T DW? - s - = dada Рә. Wi a 
> —һ CH e "1 CW AAA ds фә ӱр" y == : T = x pd 
---- LE ә UU ә ә u. 5... pes ҹ̧ : a : : : = = | : = 
` pr LLL lim. € qa - —— T E a DAL ~ E ei e a P ә, ə ә Ан нар = Nic eei 
E - - ти -—— ` ape, cm, 9, NI ~ e tn E de QN T ‘ .- ^ = e — Pam r " ” a - TR PA EE qe = ә 
лл me wa. Pape pa ə A - TER WE H - DEE E m a ə ke d. E ә ә — æ ә H oe T. E "P e пж = - 
ee KH p XA Vp s HEN ə a. yen ao oy “ T ii = — “a. “ Le ” ^ PO A < e = - - 10. 
Nee Eh ¿A “ə ah eo mp A a hn ` Ae A ад E a E x = = "R ee: Ë e ëng ә ыб -o ET ——. - P 
TT. S x 2 A d . e Geht a v pa a on = EI Sane a .. ә D ~ La "A . — pem - a P Pr * e - 
a ee pepe M Y уҹу. E v — ^ M - ^ H er A - Ko. = me T€ es ә clar «eia = 2n 9 e x 
- s. pipe pen <. < < a ew a e E TE Y dad ə... ^ bh Gab AS 7 PU d ә Le E t ml ei e = ə afp LIT d pun e дн = — pt 
= de 5 | Je T Zen KA öyə = e Dem EE fb EE tegen та KR < ^ - > Tr di A E had e vb dü * e EE ET m әәә. üsu 
2 : ES X E o — . —- va p AUN P E ти - me ә WE ~ - = w m leg Ln nai ears P ” e vio ac pr ain - ә — кии . 
^ » - d R“ qıy are Ex y 8 “A... Wa — Ech Кит a + ES e Wis T ual e “~~ Lad x” Ce oy = P * mg e мәб m sdo m PT T PP d a. -.— - жд. = e < P ~ 
= ae ae EE PN a» dS Am, par^ em a Pme m Рә ии: без" Papo em Mi de e ee oe "лн a ти. - LX E m At et WH a ән = eg e a aia sa <€ Y 
SCH Key x Ç - ə Sg S. Seem ән ҝҝ̆рә, ec Paşaya ҺР. ыр E Ly Dui ed düyə" id – y pass AA Se po = P = - - gap Fue L A EA P AO AAA ын 
ҹ = x. 9 ` ә x əə — “ e Sa - yaa p = 20. os My "P o dz “lə ep A лл ndo dl mdi ЫД. ән - 40 0-дә de бр P r T wo anm am s — ër ə c 
o. ve = a x X ee a ebe pp ya DEEN — Le < ҹа, e CA A ə md .. ara ləsi m = *– PR вы “ədə a ә. - - P LS wa 
- ^ ep EE io < ^ v ERA a o tenes ә ES s ... E ... .”c— y per aos rr — par” e e TRAD eg Pe d 
ees x pa тә, ü wen Pago qae E Lr EN "T hn te Pr? -— гр" 7 ecd» =r dür dad EE ipaa ҹә pe 
~ r Q < ^ 3 - D а da A = BAS A D p EE aman m SEET дә naga De w geg aet = 
ә ә ә ФК ә Py ә" RK deg E E “ o a, s .. ra - P Kem w ə wer t ə p. E T a m DN - nd "s d o 
| . kə = ə lege ә " P Me e di "a d P — e Pvt am e .... a ET ə мәһс дин ta — 
y - a m "e a A — a ¿A e + p pu de - BE E -— M d b a TE gë bet >si 
"— = әт" 9 N el = wd T: we = Lä RA AA AS IEA PPS As de De JO gener IO A DADA əz şı EE — qi e 
Pm tT ed eed qa -—— «x am m TT. Pv dea w = n pa v tu Ҹу eg q ao ew 
MAS ae ҹу. rera er) OS DR EEN ne - "PIS eg A | iam. gay A ada laren Ap — 
e 177 -— dame = - ës, d i E ^- AA ete ә ә ии, Ра L TO — mi ud ^ ms абы O JUN" Su a? Pa m: 
Kee WK "paya qaq - e qué ... p eos rm РӘҸРФ p әтри P" e pepper FEO RIA 
> U E al `A A alee ЛА CR oe a E ET wi p Garë ” e ..... Ps us amc m I a Pr idi .—.—.t Рә ит fÉ. s... P a ыл V anum A rte ədl 
Let EE Ehe ә = Kéi “ P Sé s ФА " po a pa Py Pg AID DS E EAN PA get tenen po 
Pa CE V biyan ^ иә ә ә ә e et Ag a — hn 4 = Ke? EE E EE E ph Date o oo ee 
< Ed ә" Aeg! d een) whet am də m dy Mu ~ ex PI ә ҺЛ P Рди E? NET E AA — e pud o A eeng ee ëch == . 
E A A de ia —— dayi, EN ame s e ET T "тут - RT TEE E E d 9 ox (bər, ei en o dinin əsə EE ADS gən s m GO qe алын s.s = 
о : KK? "PUSO SIA u (po ə wm © prat WANN ҹи, A A AS rió AO SA AR A eg 
O) i D = PA ciate teal WW Ld glo qe. x * e ef Peme pL eg L... - TT T gp >. pd pha o  — enr Ы ыы Ди 
a "pae ny mw TC a ә = ` " ә а- pee pu pups pu ya e E EH re P ə - ә ән ус ASI SN S a ae A O bua a m... = 
eO — J rap rr es Pes Pope qr ed oe Pp. e D "nr p T + PAIPA ISI Só A RR A uo agasiap ss 
ca py. 55 əvə. ılıq Re + d EI — € X e fad enfe th A gl AAA O O s a up, es af MAA A ire ARS en 
ит, Ы Zeg PL | » wa. = —. H = * KEE ae su. Sr »* ^» E > кк. "“ 2 —— s perpeti —.. PP to rs P anb 5 
о c = ` w A e D — sadi = * P — e гә ideal - PD OPA oA BP APL OLAS ә ә qiie дајә pon — ыр ДА basi 
> ANN eer ə Ҹи BE A Y iqup n A Pree pw —L n altel er e e A NO e qt ena td ap 
; "dis ya ^. pd. E EE ET = NEE a Set pou > e уҹ". — mo —-к pr A p pardos een ep ee eae agar? ORE 
= əə ə geg seg ` ^. p a Pt TTT gg P we = ҺР әј Р. ыы ыы ын = ye P TT gege, goa as, atem oum p tap dre sic ge 
O O Ve ER e A əla” "Cr əsən E rə” noui pu DY E E gegen os SS A a geg ASAS PLA pagó “— 
= qaya ən "rant sarın pp Ps. ` TS ER - w w Ру әт Ол lagna aes oe — дә PP а. a ar a UA A INR AAN AS APS SA 
O AAA rr py ad O dap dere bg nn yiq Ру" E E E ld ə E ege SARA ds apta qut ——.. PEA SE € 
sag a e TT a з ә Ы ә ды ад Qr sean qq ep tip SEPP EE RADA doo pati craint eege ee 4 SECH 
O ` APA wer rr m $n Mp Pur P “ ә. AAA A E alt ee: AS BA mp NES is SRA AR a ro A z “ə 
< P. x zo ebe: ee 9 € PC ESA — "7 P Lineas q Aba s dan Н да абад e ә poke niall DALE AEDT LOL FIO rario gi qi gu gana at. figi ad 
ah R > a, EN inis heed R m. Lud LI i Lei əl 3 Gees -== 
e | ә Pr : 
wo y Š | = 
+ . 
Mo | 
cN > 
a z 


CH, 


An 


LO 
be EE E reli di Gesi) 
parer eegen, ART äng S Gester 
Kë EE br nədi iə de? dr ap Ae 
eg EE E ee pip ды a Yw. 
pr tte 





A 
udi 


La. < 
e ët WT EE KR "T Kë Ed mé Fires em xd e 
şəh ot A A a ad RARAS “ən H e ə ye sp. <. 7 
HA DAA AAA Ҹи ти, LI 4 t A Za ° Ces? 
Ba ECK E An A cid A ə CACO an Aaa E ect tl de os 1 
- 5 == m= a "wi 





